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(57) Abstract 



A threshold estimation method in a telephony system 
estimates (SI) downlink power R* and uplink power R y . 
Echo model parameter candidates a; and bj are recursively 
determined (S2) for delayed R x values. The X>\ having the 
minimum value is selected (S3) as an echo model parameter 
bi and in the corresponding delay A (index i) is identified 
(S4) as the current echo delay. The model parameter 
candidate a* having the same delay is selected (S5) as model 
parameter a. The current echo return loss ERL is estimated 
(S6) from the estimated values of a and b. The current clip 
level is estimated (S7) from the estimated value of ERL 
and the estimated value of R, that was measured A time 
steps previously. 
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- METHOD AND APPARATUS FOR ECHO ESTIMATION AND 

SUPPRESSION 

TECHNICAL FIELD 

The present invention relates to a method and apparatus for echo estimation and 
suppression in telephony systems. 

BACKGROUND OF THE INVENTION 

Echo is a problem related to the perceived speech quality in telephony systems with long 
delays: e.g. telephony over long distances or telephony systems using long processing 
delays, like digital cellular systems The echo arises in the four-to-two wire conversion in 
the PSTN/subscriber interface (PSTN = Public Switched-Telephone Network). To 
remove this echo, echo cancellers are usually provided in transit exchanges for long 
distance traffic, and in mobile services switching centers for cellular applications. 

Due to the location of the echo canceller it is made adaptive; the same echo canceller is 
used for many different subscribers in the PSTN. This adaptation is necessary not only 
between different calls, but also during each call, due to the non-fixed nature of the 
transmission network, e.g. phase slips, three-party calls, etc. 

The main part of an echo canceller is an adaptive filter. The filter generates a replica of 
the echo, which is subtracted from the near end signal. Due to imperfect knowledge of 
the echo generating system, the estimated echo signal always contains errors. Hence, in 
practice, the echo attenuation obtained by using an adaptive filter is usually at most 
approximately 30 dB. For long time delays this attenuation is not enough, and in order to 
minimize the audible effects of these errors, a residual echo suppressor is used. The 
purpose of the echo suppressor is to further suppress the residual signal whenever this 
signal is dominated by the errors in the echo estimate. This is done by blocking the 
output of the echo canceller for certain levels of the output signal. 
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Reference [1] describes an echo canceller provided with an echo suppressor in the form 
of an adaptive center clipper. The echo estimate produced by the echo canceller is used 
to control, via signal processing means, the threshold, and thereby the clipping window, 
of this adaptive clipper. If the power of the residual signal falls below the adaptive 
threshold, the residual signal is blocked or clipped-, otherwise the residual signal is 
passed, through the adaptive clipper without modification. However, the residua |--signaT 
contains not only residual echo, but also background noise produced at the near end 
subscriber. Occasionally residual echo samples and background noise samples add 
constructively, and the resulting residual signal may therefore exceed the threshold. The 
result is undesirable sporadic transmissions of residual signals containing residual echo, 
which can be very annoying. 

Echo cancellers in transit exchanges or mobile switching centers are shared 
devices, which can be used towards any of a multitude of subscriber lines or echo 
paths. The characteristics of these echo paths can vary substantially Particularly, 
the delay, the power level and the linearity of the echo may vary This results in 
great difficulty in determining proper, fixed, values for the echo suppressor or non- 
linear processor (NLP). Hence, it is des.rable to dynamically estimate the amount of 
echo power forwarded to the echo suppressor 



Another possible cause for echo is the acoustical and mechanical coupling of the 
loudspeaker and the microphone .n telephone handset. For a mixed analog and 
digital network, such as a PSTN with analog two wire subscriber lines, this echo 
usually causes no audible disturbance. This IS so because the echo from the 
25 handset is masked by the usually much stronger electrical echo in the four-to-two 

wire conversion Furthermore, by applying an echo canceller for the latter echo, also 
the echo from the handset will be controlled 



On the other hand, in digital telephony systems even a low level of acoustical and 
mechanical cross-talk might become audible if the transmission delay is large. This 
has been recognized in the specification of several digital systems. As an example, 
the GSM specification (GSM = Global System for Mobile Communications) requires 
an echo loss of at least 46 dB between loudspeaker and microphone, as measured 
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in the Switching System. This echo loss should be measured with the handset 
earpiece sealed to the ear. and with the volume control of the telephone set to its 
maximum. However, in normal use there is usually a gap between ear and earpiece, 
which might result in a lower echo loss. Thus, although a handset may fulfill the 
GSM specification, in practice there is still a possibility that echo might be observed 
during normal use. Hence, it may be desirable to detect and cancel aslo this echo" 
from the switching system. 



SUMMARY OF THE INVENTION 

An object of the present invention is to provide an echo power estimation and 
suppression method and apparatus that avoids the above mentioned problems. 

This object is achieved by a method an apparatus in accordance with the appended 
claims. 

Briefly, the present invention is based on the finding that there is a correlation 
between the power of a delayed echo and the power of the signal generating the 
echo. This correlation is used to build a model capable of estimating echo return 
loss and delay. These parameters are used to form a clip level used as a dynamic 
threshold for controlling an echo suppressor. The method is also computationally 
efficient. 



BRIEF DESCRIPTION OF THE DRAWINGS 



The invention, together with further objects and advantages thereof, may best be 
understood by making reference to the following description taken together with the 
30 accompanying drawings, in which: 

FIG. 1 is a block diagram of an echo generating system; 

FIG. 2 is a block diagram of an echo cancellation system; 
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FIG. 3 is a block diagram of a previously known echo canceller with" srvarrable 
center clipper; . .. 

FIG. 4 is a diagram illustrating the transfer function of the center clipper of the 
echo canceller of Fig. 3; - - ; - _ 

FIG. 5 is a block diagram of an echo generating system in a cellular telephony 
system; .. . 

FIG. 6 is a scatter diagram illustrating the correlation between the power of uplink 
and downlink signals in a cellular telephony system with zero delay between the signals; 

FIG. 7 is a scatter diagram illustrating the correlation between the power of uplink 
and downlink signals in a cellular telephony system w.th delay between the signals; 

FIG. 8 is a diagram illustrating the correlation model between the power of uplink 
and downlink signals in a cellular telephony system with delay between the signals; 

FIG 9 is a block diagram of an apparatus for estimating a dynamic threshold for 
controlling an echo suppressing non-linear processor; 

FIG. 10 is a block diagram of a cellular telephony system provided with a 
threshold estimator for controlling a non-linear processor; 

FIG. 11 is a block diagram of a telephony system provided with a threshold 
estimator for controlling a non-linear processor provided after an echo canceller; 

FIG. 12 is a block diagram of a modified embodiment of a telephony system 
provided with a threshold estimator for controll.ng a non-linear processor provided after 
an echo canceller and 

FIG. 13 is a flow chan illustrating the threshold estimation method in accordance 
with the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In the following description elements performing the same or similar functions have been 
provided with the same reference designations. 

Before the invention .s described in detail, the two environments in which the invention is 
to be used will be more fully described. 
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Fig. 1 illustrates the echo generating process, in a telephony system. A subscriber A, 
called the far end subscriber below, is connected to a hybrid (a hybrid forms the 
interface between a four-wire and a two-wire connection, as is well known in the art) over 
a two-wire line. Similarly a subscriber B. called the near end subscriber below, is 
connected to another hybrid over a two-wire line. The two-wire lines transfer both 
incoming and outgoing speech signals ' Outgoing speech from far end subscriber AHs" 
transferred to near end subscriber B over the upper two-wire line in Fig. 1. Similarly 
outgoing speech from near end subscriber B is transferred to far end subscriber A on the 
lower two-wire line in Fig. 1. However, the lower two-wire line from subscriber B to 
subscriber A also contains an echo of outgoing speech from subscriber A. which the 
hybrid- at subscriber B was not able to suppress completely. Similarly the upper two-wire 
line in Fig. 1 contains echo from outgoing speech from subscriber B. 

Fig. 2 illustrates how the echo back to subscriber A is cancelled at the near end side (a 
similar arrangement is provided at the far end side). Input signal x(n), where n denotes 
discrete time, represents speech from subscriber A. The input signal x(n) is attenuated 
by the hybrid, represented by a filter 10 with transfer function H(q' 1 ), where q 1 represents 
the backward shift operator (q 1 x(n)=x(n-1)), and a summation unit 14, and the resulting 
echo signal s(n) is combined with the near end signal v(n), which may or may not contain 
near end speech, in summation unit 14. The attenuation of filter 10 is represented by the 
echo path attenuation ERL (ERL = Echo Return Loss). Thus, the resulting output signal 
y(n) contains both the near end signal and echo from the far end signal. Furthermore, 
input signal x(n) is also forwarded to an adaptive filter 12. which models the impulse 
response of the hybrid by adjusting its filter coefficients. The resulting estimate of echo 
signal s(n) is denoted S(n). This estimate is, in a summation unit 16, subtracted from 
output signal y(n), and the resulting error signal e(n) is forwarded to adaptive filter 12 for 
adjustment of the filter coefficients and to the two-wire line back to far end subscriber A. 
The coefficients of filter 12 may be adjusted in accordance with, for example, the NLMS 
algorithm (see [2]). 

Since adaptive filter 12 may not completely remove the echo, the echo canceller may be 
supplemented by an echo suppressor. Figure 3 illustrates an arrangement of this type 
described in reference (1). Here the echo suppressor has the form of a variable center 
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clipper 18/ which is controlled by the estimated echo signaf §(n): The output signal 
evcc(n) from center clipper 18 may be described-by: - ' - - 

, , f 0 : ' e(n) t = S \ s(n) \ 
I e (»> ■ ' <i(n) j > r 

where x is the variable threshold of center clipper 18 and 6 is a constant scaling factor* 
This transfer function is illustrated in Fig. 4. Thus, if s.gnal e(n) falls within" the center 
window it will be completely blocked, otherwise it will pass through center clipper 18 
unaltered. However, threshold x depends on the value of s(n). Thus, if s(ri) has a low 
magnitude, threshold x will also be low. In this case |e(n)| may actually exceed the 
threshold if signal e(n) contains significant background noise from subscriber B. Thus, 
signal e VC c(n) may contain sporadic s.gnals containing echo and background noise. The 
solut,on in accordance with the present invention described beiow will eliminate this 
problem. 

Another environment in which echo ,s generated is illustrated in fig. 5. This figure shows 
relevant parts of a digital cellular telephony system. A s.gnal x(n), called the downlink 
signal below, from a far end subscriber is forwarded to a speech coder 20. which 
transforms the speech samples x(n) in to speech code The speech code is forwarded to 
a radio unit 22 for transmission on an antenna 24 The transmitted signal ,s received on 
an antenna 26 and a radio un.t 28 of a mobile station. The received speech code is 
decoded in a speech decoder 30 and the decoded speech is forwarded to a loudspeaker 
(earpiece) 32. As indicated by the curved arrow in fig. 5, sound from loudspeaker 32 will 
also partially reach a microphone 34 of the mobile stat.on. The echo containing signal 
from microphone 34 is encoded in a speech coder 36. and reaches a speech decoder 
38 over radio unit 28. antennas 26. 24 and radio unit 22. The decoded echo s.gnal y(n), 
denoted the uplink signal below, is forwarded to an echo suppressor or NLP 40. 

The system in fig. 5 may also contain an echo canceller before echo suppressor 40. but 
since an echo canceller is a linear filter that removes only linear echo effects and speech 
coding/decoding is a strongly non-linear process, the echo canceller may also be 
omitted, as in fig. 5. 
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The problem that the present invention solves is how to control echo suppressors 18, 40 
in an efficient way. The solution will now be described in detail with reference to fig. 6- 
13. 

Fig, 6- is a scatter diagram illustrating the dependence of the power R y of the uplink 
signal. y(n) in fig. 5 omthe power R x of the downlink signal x(n) when the delay between* 
the signals has been ignored. The powers have been calculated in accordance with the 
equations: 



.V-1 



(Other power calculation methods, for example an exponential window may also be 
used.) The power is calculated every N:th speech sample (N=64, for example). The 
diagram represents many samples (480 000 speech samples or 60 seconds of speech). 
The intensity at any given point represents the number of power samples at that point. It 
is apparent from fig. 6 that there is little or no correlation between power samples if the 
delay is ignored. 

Fig. 7 is based on the same data, but in this case a delay of 208 ms has been assumed 
20 between downlink and uplink signals. The difference is evident. In this case there is a 

strong correlation between power measures R* and R y , at least for R*>-30 dBmO. For 
Rx<-30dBmO, R y is -60 dBmO regardless of the value of R*. Fig. 7 suggests that some 
kind of model may be used to describe the relationship between the uplink and downlink 
signals. As a starting point it is noted that: 



/?,.<«) = /?,(») + 

where R s represents the power of the echo and R v represents the power of the signal 
produced at the near end (mobile). By definition the power of the echo, Rs, is: 



30 
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and therefor R y may be written as: 

, >>) V ^ + /?(") 

Computing the logarithm of both sides, this may be approximated by: 

l°gk.(»))= {' og(/? - (/?)) log(/?,.(«)) > log(y? v (»))- yERL 

1(1 .-•C)log(/? r (#i))+-.C*(log(/? 1 («) - log(7?, („)) < \og(R x (n))-yERL 

where v = log(10)/10 and C = 1 is a predetermined constant The purpose of C is to 
provide a better approximation when log R v = log R, - yERL. In this approximation R v and 
ERL are considered unknown. Introduce parameters a(n) and b(n) defined by: 

a(n) = los(/<(/;)) 
b(n) = yERL + a(n) 

With these definitions the above approximation may be rewritten as: 

log(/<,.(«))= </</>) + ( ' max(0,log(7? v (//))- b(n)) 

where the parameters a(n) and b(n) have to be estimated dynamically. This model is 
illustrated graphically in fig, 8. 

The previous equation does not give explicit information regarding echo delay. To 
determine the delay, the parameters a(n) and b(n) of the model are computed for a 
number of different delays, and the delay corresponding to the smallest b(n) estimate is 
selected as an estimate of the actual echo delay. The ERL is estimated in accordance 
with ERL = (b(n)-a(n))/ Y . where a(n) and b(n) are the parameters that correspond to the 
estimated delay. The parameter estimation is preferably performed by a recursive 
algorithm. The power estimates R, and R y are computed for a rectangular window with 
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N=64, for example. Delayed estimates of the downlink power estimate R* are stored up 
to a delay corresponding- to an expected maximum echo delay: For example;- if- it known 
that the echo delay is at most 264 ms. a total of 33 delayed power estimates need to be 
stored (33x64 = 2112 samples, which corresponds to 264 ms for a sampling frequency 
of 8000 sample/s). Estimates of the parameters a and b are then computed recursively 
every Nth sample for each delay by the following LMS type algorithm: - ■ - 



e,(n) = \og(R y (n))- a, (») - C max[0, log(/?,(/i - / N)- 
o,(« + l) = a(«) + ^ <) •<?,(») 

A i (« + l) = P ,( " ) ~ //fc ' Ce ' (/,) l °8(^("-i N))>b,(n) 



otherwise 



/ = 0 • • - / 



i 0 where 

i is an index representing a delay by i N steps of discrete time, 
imax is an integer representing the maximum delay by i™ N steps of discrete time, 
a(n) and bi(n) denote the model parameter estimates corresponding to a delay of 
i samples, 

ei(n) denotes an estimation error associated with a delay of i samples, and 
\u and nb denote step length parameters of the algorithm. 

In the algorithm "log" may be implemented as the integer part of the 2-logarithm to 
further reduce complexity. 
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This algorithm is performed by block 42 in the threshold estimator 50 in fig. 9. In block 44 
the echo delay is identified as the discrete time delay that gave the smallest bi(n) value. 
This minimum value is also used as the current parameter value b(n). Thus, the delay a 
is calculated as: 

A(/7) = delay{b mn {n)) 

The determined index i and the corresponding parameter value b(n) are now used by 
block 46 to estimate the echo return loss ERL according to: 



30 
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• # 

"" ....... - y 

where a(n) .s the value ai (n) that has the same delay as b(n) (the same index i). 

Finally block 48 estimates a dynamic threshold or clip level using ; the.estimated,delay_A 
and echo return loss ERL in accordance with: 

clip level =k RA "m L A ± ■ 



where K is a constant, preferably of order 10. mtroduced to account for uncertainties in 
10 the estimate of ERL 

Preferably the functions of the blocks of threshold estimator 50 are performed by one or 
several micro/signal processors 



Fig. 10 illustrates a digital cellular telephony system similar to the system in fig. 5, but 
provided with a threshold estimator 50 in accordance with the present invention for 
controlling echo suppressor 40. As described above threshold estimator 50 estimates a 
dynamic clip level. This clip level is forwarded to echo suppressor 40. If the power of 
uplink signal y(n) falls below this clip level, the signal is assumed to be dominated by 
echo and is blocked, otherw.se it assumed to be dominated by the near end signal from 
the mobile station and is passed through echo suppressor 40 without change. If the 
signal is blocked, comfort noise may be added instead. The comparison with the 
dynamic clip level is performed on a sample by sample basis. 

The same type of echo suppressor control as described above with reference to a 
cellular telephony system may also be used in a regular PSTN. This is illustrated «n fig 
11 and 12. 

Fig. 1 1 is a block diagram of a telephony system that is provided with a threshold 
estimator 50 for controlling a non-linear processor or echo suppressor 18 provided after 
an echo canceller. In this embodiment the "downlink" signal corresponds to signal x(n) 
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from the far end subscriber, while the "uplink" signal is represented by residual echo 
signal e(n). Under these assumptions echo suppressor 18 is controlled in the same way 
as echo suppressor 40 in fig. 1 0. 

Fig. 12 is another embodiment of a telephony system that is provided with a threshold 
estimator 50 for controlling a non-linear processor or echo suppressor-18 provided after 
an- echo canceller. In this embodiment the "downlink" signal corresponds to echo 
estimate signal S( n ) instead; of far end subscriber signal x(n). The reason for this is that 
since the echo estimate has to be calculated to control the echo canceller anyway, and 
since it obviously is strongly correlated to signal x(n). it might as well be reused for the 
threshold estimation. 
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FIG. 13 is a flow chart illustrating the threshold estimation method in accordance with the 
present invention. In step S1 current estimates R, and R y are calculated. Step S2 
recursively determines model parameter candidates a and bi for delayed R x values. In 
step S3 the bi having the minimum value is selected as model parameter b. and in step 
S4 the corresponding delay a (index i) is identified as the current echo delay. The model 
parameter candidate a having the same delay is selected as model parameter a in step 
S5. Step S6 estimates the current echo return loss ERL from the estimated values of a, 
b and delay. In step S8 the current clip level is estimated from the estimated value of 
ERL and the estimated value of R, that was measured a time steps previously. Finally 
the process loops back to step S1 to determine the next estimates of R* and R y (at the 
same time the oldest value of R, is dropped). 

The present invention is a computationally efficient method for detecting and estimating 
power level and delay of a low level echo, such as the echo originating from acoustic 
cross-talk in mobile terminals or the residual echo remaining after the linear filter of an 
echo canceller. The determined parameters may be used to control the operation of a 
switch (NLP) which blocks transmission when echo is present. 

For elimination of cross-talk echo the present invention may result in great computational 
savings, since a standard echo canceller may be replaced by a dynamically controlled 
echo-suppressing switch. Due to the non-linearities of the echo path in this application, 
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the subjective performance of the proposed method is similar to the performance of a 
much more complex echo canceller. 

It will be understood by those skilled in the art that various modifications and changes 
may be made to the present invention without departure from the spirit and scope 
thereof, which is defined by the appended claims. - * 

REFERENCES 
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[2] D IM. Slock, "On the Convergence Behavior of the LMS and the Normalized LMS 
Algorithms". IEEE Transactions on Signal Processing, 41 (9):281 1-2825, September 
1993 
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CLAIMS 

- r T! Ah echo power estirriatioh method in a telephony system, characterized by 

forming an echo power estimate R y (n) of an echo signal, said echo power 
-estimate depending on an input power estimate R x (n) of an input signal in accordance 
with the approximation: 



fogfa (»))= a(//) + C max (OJogfo («))-£(«)) 

where 

10 n represents discrete time, 

a(n) and b(n) are dynamically updated model parameters, and 
C is a predetermined constant 

2. The method of claim 1, characterized by determining said model parameters a(n.) and 
1 5 b(n) in a recursive algorithm. 



3. The method of claim 2, characterized by 

a recursive algorithm in accordance with the equations: 

e, (n) = \og(R ,.(/;))- a, (/;) - C max [O, log(^ A (n - / • A') - />, (>;)] 
</,(/> + !) =«,(«) + //, e(n) 

h (/, + I ) = J h ' in) ~ C ' e - { " ] lo ^ ( " ~i-N))> b. (/,) 

otherwise 

20 

where 

N is a predetermined integer, 

i is an index representing a delay by i N steps of discrete time, 
imax is an integer representing the maximum delay by i max N steps of discrete time, 
2 5 a(n) and bi(n) denote model parameter estimates corresponding to a delay of i 

samples, 

e,(n) denotes an estimation error associated with a delay of i samples, and 
Ha and jib denote step length parameters of the algorithm; and 



/ = 0- 
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selecting the model parameter estimates a(n) and b(n) that are associated with 
the index i that gives the smallest value .ot b^as the current model parameters a(n) 
and b(n). 

4. The method of claim 3, characterized by 

calculating said input power estimate in accordance with the equation: r- -- - 

^(») = £v : ( / /-A) 

where x(n) is said input signal; and 

calculating said echo power estimate in accordance with the equation: " ' 

where y(n) is said echo signal. : 

5. The method of any of the preceding claims 1 -4, characterized by 

estimating an echo delay a as the delay of the selected current model parameter 

15 b(n). 

6. The method of any of the preceding claims 1-4. characterized by 

estimating an echo return loss ERL in accordance with the equation: 

20 ERL(») = a( " ) - h(n) 

r 

where y is log (10)/10. 

7. The method of any of the preceding claims 1-6, characterized by sa.d echo signal 
being a residual echo signal remaining after an echo canceller. 



8. The method of any of the preceding claims 1-6, characterized by said input signal 
being a downlink signal to a mobile station and said echo signal being an uplink signal 
from said mobile station. 
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9. Ahfecho suppression method in a telephony system; characterized 'by: 

— ■ —forming ah echo power estimate" R/n) of an echo' signal,' said echo" power 
estimate depending on an input power estimate R,(n) of an input signal in accordance 
with the approximation: 

•log(/?..("))=a(») + C max(0.log(/? v (n))-6(/0) '■ . - V ." - 

where n represents discrete time, a(n) and b(n) are dynamically updated model 
parameters, and C is a predetermined constant; 

estimating an echo delay a by the model parameter b(n); 

estimating an echo return loss ERL in accordance with the equation: 

IZRI = ai ") 

Y 

where y is log (10)/10; 

estimating a dynamic clip level in accordance with the equation: 

clip level = (/? ~ A) 
10'" 

where k is a predetermined constant, and; 

suppressing said echo signal whenever its power falls below said clip level 

10. The method of claim 9, characterized by determining said model parameters a(n) 
and b(n) in a recursive algorithm. 

1 1 The method of claim 10, characterized by 

a recursive algorithm in accordance with the equations: 

e, («) = log(rt»)- a, (/») - r max[0, logfo(« - i ■ N)- b,{n)] 
+ a, + . , = 0 ..., 

b, (n + 1) = {*• (,?) " * ' C • e ' (n) ,0 b(*.<« ~i ■"))>*>, (») 
\b,( n ) otherwise 
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where 

N is a predetermined integer, 

i is an index representing a delay by i N steps of discrete time, 
w is an integer representing the maximum delay by i™ N steps of discrete time, 
a(n) and b f (n) denote model parameter estimates corresponding to a delay of 
samples. 

e(n) denotes an estimation error associated with a delay of i samples, and 
M* and ^ denote step length parameters of the algorithm; and z 
selecting the model parameter estimates a(n) and b^DMhat^re associated with 

the index i that gives the smallest value of b i( n) as the current model parameters a(n) 

and b(n). 

12. The method of claim 1 1 . characterized by 

calculating said input power est.mate in accordance with the equation: 

V-l 

where x(n) is said input signal; and 

calculating said echo power estimate in accordance with the equation: 

.V I 

where y(n) is said echo signal. 

13. The method of any of the preceding claims 9-11, characterized by said echo signal 
being a residual echo signal remaining after an echo canceller. 

14. The method of any of the preceding claims 9-11, characterized by said input signal 
being a downlink signal to a mobile station and said echo signal being an uplink signal 
from said mobile station. 

15. An echo power estimation apparatus in a telephony system, characterized by 

means (42) for forming an echo power estimate R y (n) of an echo signal, said echo 
power estimate depending on an input power estimate R,(n) of an input signal in 
accordance with the approximation: 
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log(/? v («))= a(n) + C max(0, log(* v (/,)) - b(n)) 

where 

n represents discrete time, 

a(n) and b(n) are dynamically updated model parameters, and 
C is a predetermined constant. 

16. The apparatus of claim 1 5, characterized by means (42) for recursively determining 
said model parameters a(n) and b(n). 

1 7. The apparatus of claim 16, characterized by 

means (42) for- recursively determining said model parameters a(n) and b(n) in 
accordance with the equations: 

e,(n) = log(R y (n))- a (») - C max[o, log(/?> - i N)- b{n)) 
d,(n + I) = a i («) + // o e,(n) 



/>,(" + !) = \ h,{n) ~ - U o~ C e ' (/,) lo g(^<" ~i N))>b, (/;) 
lA(«) otherwise 



/ = 0 — / 

max 



where 

N is a predetermined integer, 

i is an index representing a delay by i-N steps of discrete time, 
imax is an integer representing the maximum delay by \ max -N steps of discrete time, 
ai(n) and bi(n) denote model parameter estimates corresponding to a delay of i 
samples, 

ei(n) denotes an estimation error associated with a delay of i samples, and 

M* and ub denote step length parameters of the algorithm; and 

means (42) for selecting the model parameter estimates a(n) and b s (n) that are 
associated with the index i that gives the smallest value of bi(n) as the current model 
parameters a(n) and b(n). 

1 8. The apparatus of claim 1 7, characterized by 

means (42) for calculating said input power estimate in accordance with the 
equation: 
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v-i 



where x(n) is said input signal: and 

means (42) for calculating said echo DO w*r * 
equation: 6Stimate ,n accor dance with the 

V -I 

*,.(«) = 2> : <»-*) 

where y(n) is said echo signal. 

1 9. The apparatus of any of the preceding claims 1 5-1 8. characterized by 

-el P ™i: 6St,matin9 ^ 6Ch0 - * - - ^ - - — -nt 

20 The apparatus of any of the preced.ng c.a,ms 1 5-18. characterized by 

means (46) for estimating an echo return «= D1 
equation: SS ERL ,n ac «^dance with the 



ERUm = OP)'/*") 
Y 

where y is log (10)/10. 

21 An echo suppressor in a telephony system, charactenzed by 

means (42) for forming an echo power estimate R nf 
Power est,mate depend.ng on an input pcT^ p 

accordancew,ththeapprox,mat,on: * ^ ° f ^ «** * 

"osK^^af/O + r-maxCo.log^^))-^)) 

where n represents discrete time afn) and h#«\ 

« uine, a(n) and b(n) are dynamically undated mnrioi 

parameters, and C is a predetermined constant 

means (44) for estimating an echo de.ay A by the mode, parameter b(n) 
means (46) for estimating an echo return loss ERL in ^ H 

equation: accordance with the 
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ERL{n) = a <">-*<"> 

where y is lag ( 1 0)/1 0; 

means (48) for estimating a dynamic clip level in accordance with the equation: 



clip level(n) =y ^ (/7 ^ A) 
10 ,0 . 

where k is a predetermined constant, and; 

a non-linear processor (18) for suppressing said echo signal whenever its power 
falls below said clip level. 

22. The echo suppressor of claim 21. characterized by means (42) for recursively 
determining said model parameters a(n) and b(n). 

23. The echo suppressor of claim 22. characterized by 

means (42) for recursively determining said model parameters a(n) and b(n) in 
accordance with the equations: 

! 

e\n) = log(/? y (/i))- a t {n) - C max[o, \og{R t {n - i N)- b t (n)] 

a(/7 + 1) = a,(») + ^ * (») \ i =0'-' m >s 

, . lv \b(n)-u h C e.(n) \oz(R x (n -/ - N))> b(n) I 

{b t (n) oihemise ) 



where 

N is a predetermined integer, 

i is an index representing a delay by i N steps of discrete time, 
imax is an integer representing the maximum delay by ima* N steps of discrete time,— 
ai(n) and bi(n) denote model parameter estimates corresponding to a delay of 
samples, 

e;(n) denotes an estimation error associated with a delay of i samples, and 
V* and denote step length parameters of the algorithm; and 
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means (42) for sereBTing the model parameter estimate^k(n) and b*(n) that are 
associated with the index i that gives the smallest value of b*(n) as the current model 
parameters a(n) and b(n). 

24. The echo suppressor of claim 23, characterized by 

means (42) for calculating said input power estimate in accordance with the 
equation: 

/? x («) = 2y(«-*) 

where x(n) is said input signal: and 

means (42) for calculating said echo power estimate in accordance with the 
equation: 

/v in) = V y : {n-k) 

where y(n) is said echo signal 
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